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Abstract 

Background: Triacylglycerides (TAGs) are a class of neutral lipids that represent the most important storage form of 
energy for eukaryotic cells. DGAT (acyl-CoA: diacylglycerol acyltransferase; EC 2.3.1.20) is a transmembrane enzyme 
that acts in the final and committed step of TAG synthesis, and it has been proposed to be the rate-limiting enzyme 
in plant storage lipid accumulation. In fact, two different enzymes identified in several eukaryotic species, DGAT1 and 
DGAT2, are the main enzymes responsible for TAG synthesis. These enzymes do not share high DNA or protein 
sequence similarities, and it has been suggested that they play non-redundant roles in different tissues and in some 
species in TAG synthesis. Despite a number of previous studies on the DGAT1 and DGAT2 genes, which have 
emphasized their importance as potential obesity treatment targets to increase triacylglycerol accumulation, little is 
known about their evolutionary timeline in eukaryotes. The goal of this study was to examine the evolutionary 
relationship of the DGAT1 and DGAT2 genes across eukaryotic organisms in order to infer their origin. 

Results: We have conducted a broad survey of fully sequenced genomes, including representatives of Amoebozoa, 
yeasts, fungi, algae, musses, plants, vertebrate and invertebrate species, for the presence of DGAT1 and DGAT2 
gene homologs. We found that the DGAT1 and DGAT2 genes are nearly ubiquitous in eukaryotes and are readily 
identifiable in all the major eukaryotic groups and genomes examined. Phylogenetic analyses of the DGAT1 and 
DGAT2 amino acid sequences revealed evolutionary partitioning of the DGAT protein family into two major DGAT1 
and DGAT2 clades. Protein secondary structure and hydrophobic-transmembrane analysis also showed differences 
between these enzymes. The analysis also revealed that the MGAT2 and AWAT genes may have arisen from 
DGAT2 duplication events. 

Conclusions: In this study, we identified several DGAT1 and DGAT2 homologs in eukaryote taxa. Overall, the data 
show that DGAT1 and DGAT2 are present in most eukaryotic organisms and belong to two different gene families. 
The phylogenetic and evolutionary analyses revealed that DGAT1 and DGAT2 evolved separately, with functional 
convergence, despite their wide molecular and structural divergence. 



Background 

Triacylglycerols (TAGs), fatty acyl ester derivatives of 
glycerol, are a class of neutral lipids that represent the 
most important storage form of energy for eukaryotic 
cells [1,2]. In a number of plant species, TAGs are 
major storage lipids that accumulate in developing 
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seeds, petals, pollen grains, and fruits [3]. Plant oils have 
been used for human consumption and have become 
important renewable resources as biofuels [4,5]. The 
enzymatic machinery for the formation of TAGs is 
located in the endoplasmic reticulum (ER). TAGs can 
then accumulate as oil droplets in the cytoplasm or in 
specialized oil storage bodies [6], which are generated 
through budding of the outer ER membrane [7]. A sub- 
stantial part of TAG synthesis is performed by enzymes 
of the Kennedy pathway, which sequentially transfer 
acyl chains from acyl-CoAs to sn-1, -2 and -3 positions 
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of glycerol backbones [8,9]. DGAT (acyl-CoA: diacylgly- 
cerol acyltransferase; EC 2.3.1.20) is a transmembrane 
enzyme that functions in the final step of TAG bio- 
synthesis, catalyzing the acylation of sn-l,2-diacylgly- 
cerol (DAG) at the sn-3 position using an acyl-CoA 
substrate. DGAT has been proposed to be the rate-lim- 
iting enzyme in plant storage lipid accumulation [10,11]. 
Consequently, DGAT is considered a key enzyme for 
biotechnological purposes; it might be utilized to 
increase oil content in oleaginous plant species [12-15]. 
DGAT overexpression causes a net increase in seed oil 
content in Arab idop sis [16]. Additional evidence is also 
available from studies of Glycine max [17], Brassica 
napus [18] and Zea mays [19]. TAG can also be formed 
by an acyl-CoA-independent pathway. In this pathway, 
the phospholipid:diacylglycerol acyltransferase (PDAT; 
EC 2.3.1.158) enzyme responsible for TAG synthesis, 
using phosphatidylcholine (PC) as the acyl donor, in 
which the transfer of an acyl group from the sn-2 posi- 
tion of PC to the sn-3 position of DAG yields TAG and 
sn-1 lyso-PC [20-22]. 

Different types of DGAT enzymes have been identified 
in several species. DGAT1 and DGAT2 are two of the 
enzymes that are responsible for the main part of TAG 
synthesis in most organisms, and they have been studied 
in many eukaryotic organisms. Another group of DGAT 
enzymes is represented by the bifunctional DGAT/wax 
ester synthase (ADP1) from Acinetobacter calcoaceticus 
[23]. Homologs of ADP1 have also been characterized 
in Petunia [24] and Arabidopsis [25]. A soluble form of 
DGAT has been identified in peanut cotyledons [26], 
and in Euonymus alatus, a distinct type of DGAT was 
recently identified that is responsible for the synthesis of 
3-acetyl-l,2-diacyl-sn-glycerols (acTAG), an unusual 
triacylglycerol [27]. The first cloned DGAT gene was 
isolated from mice; it corresponds to a member of the 
DGAT1 family, which has high sequence similarity with 
the sterol:acyl-CoA acyltransferase (ACAT; EC: 2.3.1.26). 
DGAT1 is structurally related to the ACATs, with the 
divergence in its amino acid sequence conferring its 
substrate specificity to DAG. Both enzymes belong to a 
large family of membrane-bound O-acyltransferases 
(MBOAT) [28,29]. The DGAT2 gene was initially identi- 
fied in the oleaginous fungus Morteriella ramanniana 
[30]. In addition to DGAT2, this family includes acyl- 
CoA:monoacylglycerol acyltransferase- 1 EC:2.3.1.22 
(MGAT1) [31], MGAT2 [32,33], MGAT3 [34], and acyl- 
CoA wax-alcohol acyltransferase (AWAT, EC:2.3.1.75) 
[35,36]. Interestingly, DGAT2 enzyme family members 
do not show DNA or protein sequence similarities with 
DGAT1. Furthermore, DGAT1 proteins are larger than 
DGAT2 and possess 6 to 9 transmembrane domains, 
while DGAT2 has only one or two [2]. Thus, it has 
been suggested that DGAT1 and DGAT2 play non- 



redundant roles in different tissues and species in TAG 
synthesis [37-40]. Homologs of the DGAT1 and DGAT2 
genes have been identified in several eukaryotic organ- 
isms [15,30,36,41-49]. A study of comparative genomics 
and proteomics in vertebrates suggested that the 
DGAT1 and DGAT2 gene families have evolved sepa- 
rately during vertebrate evolution [50]. Despite many 
studies of the DGAT1 and DGAT2 gene families, which 
emphasize the importance of these enzymes in biotech- 
nology for increasing the accumulation of triacylglycer- 
ols [16,51,52], little is known about their evolutionary 
origins in eukaryotes apart from vertebrates. 

The increasing availability of whole genome sequences 
and the annotation of genes from a wide range of phyla 
enable deeper phylogenetic analysis of gene families to 
provide important contextual insight into their present 
day diversity and evolution. The goal of this study was to 
examine the evolutionary relationship of the DGAT1 and 
DGAT2 genes across eukaryotic organisms in order to 
determine their origin and to verify if these genes evolved 
separately in eukaryotic organisms. We performed a 
comparative analysis integrating phylogenetic, computa- 
tional and structural approaches to the DGAT1 and 
DGAT2 genes, including analysis of their relationship 
with the ACAT, MGAT and AWAT genes in a large 
number of fully sequenced eukaryotic genomes, repre- 
senting a broad array of taxonomic groups. Most gen- 
omes contain both the DGAT1 and DGAT2 genes, 
though DGAT1 is absent in yeast and basidiomycetes. 
Our findings suggest that these two enzymes evolved 
separately since the emergence of eukaryotes and may 
have different origins. 

Methods 

Data sources and sequence retrieving 

DGAT1 and DGAT2 gene and protein sequences were 
obtained by performing blast searches (blastp, blastx and 
tblastx) of the Protein and Genome databases with the 
default parameters and an e- value threshold of 1.0 E-50 at 
the NCBI (National Center for Biotechnology Informa- 
tion-http://www.ncbi. nlm.nih.gov) and the completed 
genome projects database at the JGI (Joint Genome Insti- 
tute-http://www.jgi.doe.gov). The cDNA sequences encod- 
ing DGAT1 and DGAT2 from Arabidopsis (AT2G19450 
and AT3G51520) were used as queries in the blastx and 
tblastx programs to search for DGAT1 and DGAT2 from 
different plant species. The predicted peptides were 
further used in blastp analysis to double check their pre- 
dicted identities. The mouse DGAT1-2 (NM_010046 and 
NM_026384), ACAT1-2 (NP_033256 and NP_666176), 
MGAT1-2 (NM_026713 and NM_177448) and AWAT1-2 
(NM_001081136 and NM_177746) cDNA sequences were 
used as queries in the blastx and blastp to search for 
homologous sequences. The genome databases of yeast 
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http://www.yeastgenome.org/, fungi and insects http:// 
genome.jgi-psf.org/ were also analyzed using a methodol- 
ogy similar to that described for the plant and vertebrate 
sequences. Sequences from representative eukaryotic spe- 
cies belonging to plant monocotyledons (Oryza sativa, 
Sorgum bicolor, Br achy podium distachyon), eudicots 
(Arabidopsis thaliana, A. lyrata, Brassica napus, Ricinus 
comunis, Manihot esculenta, Populus trichocarpa, Prunus 
persica, Medicago truncatula, Glycine max, Vitis vinifera, 
Cucumis sativus, Mimulus guttatus), mosses (Selaginella 
moellendorfii, Physcomitrella patens), microalgae (symbio- 
tic Chlorella sp., free-living Coccomixa sp.), mammals 
{Homo sapiens, Pan troglodytes, Bos taurus, Mus musculus, 
Rattus norvegicus), bird (Gallus gallus), amphibian {Xeno- 
pus tropicalis), fish (Danio rerio), invertebrates {Apis melli- 
fera, Drosophila melanogaster), yeast {Saccharomyces 
cereviseae, Candida albicans), fungi {Aspergillus niger, 
Alternaria brassicicola, Laccaria bicolor, Schizophylum 
commune, Agaricus bisporus) and Amoebozoa {Dictyoste- 
lium discoideum) were selected. Additional File 1 provides 
a detailed description of the proteins used and the corre- 
sponding accession numbers. Taxa terminologies are 
abbreviated using the first letter of the genus and two 
letters of the species name (e.g., Ath corresponds to Arabi- 
dopsis thaliana). 

Sequence Alignments 

Nucleotide and protein sequences were aligned using the 
Alignment Explorer/CLUSTALW [53] implemented in 
Molecular Evolutionary Genetics Analysis (MEGA ver- 
sion 4.0) [54]. The multiple alignments were manually 
inspected and edited, and only positions unambiguously 
aligned were included in the final analysis. The final data- 
set included a total of 151 sequences from 43 species and 
344 amino acid positions. Because the DGAT2 gene was 
first identified in the oleaginous fungus Morteriella 
ramanniana [30], we also included these sequences in 
the analysis. In Tropaeolun majus, DGAT1 seems to be 
the sole acyl-CoA-dependent DGAT [15], and this 
sequence was included in the analysis. We also included 
in this dataset sequences of other types of DGAT 
enzymes to compare with DGAT1 and DGAT2, such as 
the WSD protein sequence from Arabidopsis [25] and 
Acinetobacter sp. [55], soluble DGAT from Arachis hypo- 
gaea [26] and a distinct DGAT (DAcT) from Euonimus 
alatus [27]. 

Phylogenetic and structural analyses 

The phylogenetic analysis was reconstructed after pro- 
tein sequence alignments using three different and inde- 
pendent approaches: the neighbor-joining (NJ), the 
Bayesian and the Maximum-likelihood (ML) methods. 
The NJ method was performed with MEGA 4.0 [54]. 
The molecular distances of the aligned sequences were 



calculated according to the poisson correction model. 
All gap and missing data in the alignments were 
accounted for by pairwise deletion. Branch points were 
tested for significance by bootstrapping with 1000 repli- 
cations. The ML analysis was performed using 
PhyML3.0 [56]. Branch points were tested for signifi- 
cance by bootstrapping with 100 replicates. Bayesian 
analysis was conducted in MrBayes 3.1.2 [57] with the 
mixed amino acid substitution model plus gamma and 
invariant sites. Two independent runs of 10,000,000 
generations each with two Metropolis-coupled Monte 
Carlo Markov chains (MCMCMC) were run in parallel 
(starting each from a random tree). Markov chains were 
sampled every 100 generations, and the first 25% of the 
trees were discarded as burn-in. The remaining ones 
were used to compute the majority rule consensus tree 
(MrBayes command allcompat), the posterior probability 
of clades and branch lengths. Convergence of the two 
runs was assessed by checking the average standard 
deviation of split frequencies (below 0.01) and the 
Potential Scale Reduction Factor (PSRF, very close to 
1.00 for all parameters). 

The structural organization of the DGAT genes was 
determined after alignment of genomic DNA and cDNA 
and EST sequences. Genomic sequences were also ana- 
lyzed in the FGENESH gene structure prediction program 
http://www.softberry.com/t58]. Predicted transmembrane 
structures were obtained using the transmembrane predic- 
tion server TMHMM-2.0 http://www.cbs.dtu.dk/services/ 
and SMART database http://smart.embl-heidelberg.de/ 
with the complete protein sequences. 

Divergence estimates 

Pairwise synonymous (Ks) and non-synonymous {Ka) 
numbers of substitutions corrected for multiple hits 
were calculated using the DnaSP software (DNA poly- 
morphism analysis) [59]. 

Cloning and Yeast Complementation 

Leaves from Castor bean {Ricinus communis cv.EBDA- 
MPB1) two month-old plants were used for total RNA 
extraction using the Nucleospin RNA Plant kit (Macherey- 
Nagel) following the manufacturer's instructions. First 
strand cDNA synthesis was performed with oligo(dT) pri- 
mers and M-MLV reverse transcriptase (Promega). The 
complete open reading frames from castor RcDGATl 
(AY366496) and RcDGAT2 (EU391592) genes were ampli- 
fied using Pfx50 polymerase (Invitrogen) with the following 
primers: 5'-CACC TCTAGA ATGACGATTCTCGAAAC 
GCC-3' and 5'-CTCGAGTCAGTTCCCATCGCGATT- 
CATT-3 for RcDGATl and 5'-CACC TCTAGA ATGGGG- 
GAAGAAGCGAATCA-3' and 5'- CTCGAG TCAAAGA 
ATTTCAAGTGTAAGGTCTGC-3' for RcDGAT2. Xbal 
and Xhol sites added for subsequent manipulations are 
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underlined. The resulting blunt ended PCR fragments were 
cloned into pENTR/D-TOPO (Invtrogen) following the 
manufacturers instructions. For complementation experi- 
ments in yeast, the cloned cDNA fragments were trans- 
ferred to pVT103-U [53] as Xbal/Xhol inserts, and the 
following resulting plasmids were introduced into yeast 
(Saccharomyces cerevisiae) strains [60]: arel are2 mutant 
(HI 112), and a quadruple mutant strain arel are2 Irol 
dgal (H1246, kindly donated by Dr. Sten Stymne, Swedish 
University of Agricultural Sciences). 

Yeast growth and lipid analysis 

Single colonies of vector-only, RcDGATl or RcDGAT2 
transformed yeast in the HI 112 and HI 246 mutant back- 
grounds were inoculated in 10 mL of synthetic complete 
medium with Ura omitted and grown for 72 h at 28°C. 
Cultures were centrifuged for 4 min at 10,000 g. Cell pel- 
lets were twice extracted with 1.0 mL of chloroform- 
methanol solution (10:1, v/v) after intense agitation for 
2 min with 0.1 mm glass-beads. The chloroform phase 
was collected and evaporated under N 2 flow at room tem- 
perature. Lipids were analyzed by thin-layer chromatogra- 
phy (TLC) on silica-gel 60 (TLC Aluminium Sheets, 
Merck, Darmstadt) using the hexane/diethyl ether/acetic 
acid (90:10:1 v/v) solvent system, which is specific for 
resolving neutral lipids [61]. Lipids were revealed after 
staining with a 0.03% solution of Coomassie Brilliant Blue 
R250 (Sigma) in 20% methanol [62]. Triolein, cholesterol 
and cholesterol-ester reference lipids, used for TLC, were 
all obtained from Sigma. 

Results 

Identification of DGAT homologous sequences 

We have conducted a broad survey of fully sequenced gen- 
omes, including representatives of Amoebozoa, yeasts, 
fungi algae, mosses, plants (monocots and eudicots), verte- 
brate (mammalian, fish and amphibian) and invertebrate 
(insect), for the presence of DGAT1 and DGAT2 gene 
homologs. The full list of genes resulting from this analysis 
is summarized in Additional File 1. A set of 36 completely 
sequenced and 7 non-sequenced genomes were investi- 
gated. A total of 151 sequences were identified, of which 
91 were DGATs (DGAT1 and DGAT2), 56 corresponded 
to ACAT, MGAT and AWAT and the remaining were 
sequences of other types of DGATs. These other DGAT 
types include a soluble DGAT from peanut, a distinct 
DGAT (DAcT) from Euonymus alatus and a bifunctional 
DGAT/wax synthase (WSD) from A. thaliana and Acine- 
tobacter sp. 

It was verified that the DGAT1 and DGAT2 genes are 
almost ubiquitously found in eukaryotes, being readily 
identifiable in all the major eukaryotic groups and in all 
genomes examined with the notable exception of yeast 
(S. cereviseae and Candida albicans) and Basidiomycetes 



fungi (Laccaria bicolor, Schizophillum commune and 
Agaricus bisporus), which do not contain a DGAT1 
gene. Most organisms possess a single copy of each 
DGAT gene. The remarkable exception to this general 
pattern occurs in some members of green plants, reveal- 
ing the possible existence of lineage-specific duplica- 
tions. Homologs to ACAT genes were found in yeast, 
fungi, vertebrates (mammalian, fish and amphibian) and 
invertebrates (insects) but not in plants. The MGAT 
and AWAT genes were found only in vertebrate taxa. In 
Amoebozoa species (Dictyostelium discoideum), four 
putative DGAT sequences were identified. 

Evolutionary relationship of the DGAT gene family in 
Eukaryotes 

The 151 eukaryote DGAT amino acid sequences and 
related ACAT, MGAT and AWAT enzymes were used 
for the reconstruction of phylogenetic trees. A total of 
344 positions were included in the final dataset. Phylo- 
genetic analyses of the DGAT1 and DGAT2 amino acid 
sequences resulted in a well-resolved tree. Results from 
Bayesian NJ and ML tree analysis produced similar 
topologies and revealed evolutionary partitioning of the 
DGAT protein family into two major clades: DGAT1 
and DGAT2 (Figure 1 and Additional Files 2 and 3). 
Vertebrate and invertebrate ACAT1-2 and yeast and 
fungi ARE proteins grouped within the DGAT1 clade 
(Figure 1 and Additional Files 2 and 3). The MGAT and 
AWAT protein grouped into the DGAT2 clade. 
Although the ACAT and DGAT1 proteins are closely 
related in phylogenetic trees, the ACAT sequences from 
yeast, fungi, vertebrates and invertebrates clustered sepa- 
rately from DGAT1 sequences. Interestingly, DGAT2 
sequences from vertebrates do not form a monophyletic 
group with the DGAT2 homolog sequences from yeast, 
fungi, invertebrates and plants. Instead, the DGAT2 
sequences from vertebrates appear more closely related 
to the vertebrate MGAT2 and AWAT clades, suggesting 
that these genes may have arisen from DGAT2 duplica- 
tion events. Inside both DGAT1 and DGAT2 clades 
from plants, monocotyledons and eudicots grouped 
separately. Two of the four DGATs sequence homologs 
found in Amoebozoa (denominated A and B) grouped 
into the DGAT1 and DGAT2 clade, respectively, sug- 
gesting that DGAT1 and DGAT2 diverged very early in 
eukaryotic evolution. EaDAcT (E. alatus), AtSAT (A. 
thaliana), a soluble DGAT from peanut and WSD from 
A. thaliana and Acinetobacter sp. grouped separately 
from DGAT1 and DGAT2, but with low bootstrap and 
posteriori probability values. This result indicates that 
these types of DGAT enzymes are phylogenetically 
divergent from DGAT1 and DGAT2. Genes encoding 
EaDAcT, SAT, WSD and soluble DGAT enzymes have 
not yet been identified in other organisms. 
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Figure 1 Phylogenetic relationships between DGAT1 and DGAT2 gene families reconstructed by the Bayesian method. A total of 151 
protein sequences from eukaryotic species and 344 sites were included in the analyses. The posteriori probabilities are labeled above the 
branches. Only values higher than 0.5 are presented. Numbers within brackets correspond to the number of species within each group. The 
enzymes are represented by different colors in the phylogenetic tree: DGAT1 (red), DGAT2 (green), ACATs (purple), MGAT (dark blue) and AWAT 
(light blue). Toxo terminologies are abbreviated using the first letter of the genus and two letters of the species name: Alternaria brassicicolo 
(Abr), Arabidopsis tholiona (Ath), Dictyostelium discoideum (Ddi), Euonymus alotus (Eal), Homo sapiens (Hsa), Mortierella ramanniana (Mra), 
Thalassiosiro pseudonana (Tps), Xenopus tropicolis (Xtr). 



Structural organization and transmembrane domain 
prediction 

To further investigate the relationships between the 
DGAT1 and DGAT2 gene families, we compared the 
structural organization of the DGAT genes, when gene 
and corresponding cDNA were available. Comparisons 
were performed between DGAT1 and DGAT2 genes 
from plants and vertebrates (mouse and human). A 
detailed comparison of DGAT1 and DGAT2 genes across 
plant species and vertebrates (Figure 2) revealed a high 
degree of conservation in gene structure within each type 
of DGAT (DGAT1 and DGAT2). However, the structure 
of the genes encoding the DGAT1 protein is distinct 
from that of the DGAT2 gene. The DGAT1 genes harbor 
15 to 17 exons, while DGAT2 genes have 8 to 9 exons. In 
higher plants, DGAT1 proteins are encoded by genes that 
are generally comprised of 16 exons and 15 introns 
(Figure 2). An exception to this pattern was observed for 
M. esculenta, which probably lost its 7th exon and has 



15 exons. The DGAT2 genes are generally comprised of 
9 exons, with the exception of Arabidopsis, which has 8. 
The vertebrate ACAT has 15 exons, while MGAT has 8 
exons. Another important finding was that the structures 
of the DGAT1 and ACAT genes are more related to each 
other than to the structure of genes encoding DGAT2, 
which is more related to the MGAT gene (Figure 2). 

Predictions of transmembrane (TrM) structures were 
performed, and comparisons were made between plant 
DGAT1 and DGAT2 and mouse DGAT1, DGAT2, 
ACAT, MGAT and AWAT (Figure 3 and 4). For each 
plant DGAT1 subunit examined, 8 and 9 regions con- 
taining highly probable TrM sequences were predicted 
for monocotyledons and eudicots, respectively. For 
DGAT2, only 2 regions containing highly probable TrM 
sequences were detected in both monocotyledons and 
eudicots (Figure 3). The comparisons among mouse 
DGAT1, DGAT2, ACAT, MGAT and AWAT subunits 
revealed that DGAT1 is more closely related to ACAT, 
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Figure 2 Structural organization of plant DGAT1 and DGAT2 genes and vertebrate ACAT\ MGAT and AWAT. Exon sequences are 
represented as simplified boxes. The size of each exon (in bp) is given. Bars represent introns. Black boxes represent DGAT1 and ACAT exons 
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which has 8 high probable TrM regions, while DGAT2 is 
more closely related to the MGAT and AWAT proteins, 
with only 1 TrM region. Interestingly, S. cerevisae acyl- 
CoA sterol acyltransferase (coded by the ARE gene) 
includes 9 TrM regions, as found in plant DGAT1 and in 
mouse DGAT1 and ACAT. Further, S. cereviceae Dgal 
(homologue to DGAT2) possesses only one TrM region, 
as found in plant DGAT2 and mouse DGAT2, MGAT 
and AWAT (Figure 4). Moreover, both ACAT and 
DGAT1 proteins have an MBOAT (membrane bound O- 
acyltransferase) domain (Figure 3 and 4), while DGAT2, 
MGAT and AWAT do not. These results support the 
hypothesis that DGAT1 and DGAT2 belong to different 
gene families and evolved separately during eukaryote 
evolution, as is demonstrated by the phylogenetic tree. 
Additional Files 4 and 5 show the multiple sequence 
alignment of DGAT1 and DGAT2 protein sequences, 
respectively. In these alignments, we observed similarities 



between DAGT1 and ACAT and DGAT2 with MGAT 
and AWAT protein sequences. 

Divergence estimates 

A comparison of the rate of non-synonymous substitu- 
tions (K a ) to the rate of synonymous substitutions (K s ) 
was performed; such comparison can be used as an indi- 
cator of selective pressure acting on a protein-coding 
gene [63,64]. We calculated KJ 7<T 5 ratios for vertebrate 
DGAT1, DGAT2, ACAT, MGAT and AWAT and for 
plant DGAT1 and DGAT2. From these analyses it can 
be noted that distinct selective pressures have been 
operating on the DGAT1 and DGAT2 genes, as differ- 
ent rates were obtained (Additional file 6 and 7). Com- 
parisons between vertebrate DGAT1 and ACAT showed 
high /<r a //<T 5 ratios, while low /<r a //<r 5 ratios were observed 
between DGAT2 with MGAT and AWAT from verte- 
brates (Additional file 6). The comparisons between 
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Figure 3 Predicted transmembrane domain for plants A. thaliana (Ath), B. napus (Bna), G. max (Gma), M. truncatula (Mtr), R. comunis 
(Rco), P. tnchocarpa (Ptr), 0. sativa (Osa), and Z. mays (Zma) DGAT1 and DGAT2 sequences. See Additional File 1 for sources. The 
TMHMM web tools of the Center for Biological Sequence Analysis, Technical University of Denmark TMHMM Server plots the probability of the 
ALDH sequence forming a transmembrane helix (0-1.0 on the y-axis) (shown in red for the relevant amino acid sequences) http://www.cbs.dtu. 
dk/services/TMHMM/. Eight and nine predicted transmembrane helices are identified for DGAT1 sequences, while two to four transmembrane 
helix (or helices) were observed for DGAT2 sequences; regions of all sequences predicted to be located inside or outside the membrane are 
shown in blue and pink, respectively. The Pfam:Mboat is represented by a dotted line to DGAT1, ACAT1 and ACAT2. 
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Figure 4 Predicted transmembrane domain for Mouse (Mmu) DGAT1, DGAT2, ACAT1, ACAT2, MGAT1, and AWAT1 sequences. See 

Additional File 1 for sources. The TMHMM web tools of the Center for Biological Sequence Analysis, Technical University of Denmark TMHMM 
Server plots the probability of the ALDH sequence forming a transmembrane helix (0-1.0 on the y-axis) (shown in red for the relevant amino 
acid sequences) http://www.cbs.dtu.dk/services^MHMM/. Eight predicted transmembrane helices are identified for each of the DGAT1, ACAT1 
and ACAT2 sequences; in contrast, only one major transmembrane helix (or helices) was observed for DGAT2, MGAT1, MGAT2, AWAT1 and 
AWAT2 sequences; regions of all sequences predicted to be located inside or outside the membrane are shown in blue and pink, respectively. 
The Pfam:Mboat is represented by a dotted line to DGAT1, ACAT1 and ACAT2. 



DGAT1 and DGAT2 in both vertebrates and plants 
revealed a high level of KJK s ratios (Additional file 6 
and 7). We also compared DGAT1 and DGAT2 
between eudicots and monocots and between oilseed 
and no oilseed plants and verified low K a /K s ratios, indi- 
cating that they are under purifying selection. 

DGAT1 and DGAT2 and triacylglycerol biosynthesis 

Two different yeast mutants that are defective in TAG 
biosynthesis were used to determine if both castor bean 
DGAT orthologs are able to complement the mutated 
enzymes and allow the synthesis and accumulation of 
TAG (Figure 5). TAG synthesis was undetectable in both 



mutant strains HI 112 (arel are2) and H1246 (arel are2 
Irol dgal) carrying the empty expression vector. Only 
DGAT1 was able to reestablish TAG synthesis in the 
HI 112 strain and H1246 yeast mutants. The equivalence 
of lipid amounts among samples was confirmed by the 
homogeneous amount of mono-alkyl-diacyl-glycerol. 

Discussion 

It has long been understood that DGAT enzymes play 
important roles in triacylglycerol biosynthesis. DGAT is 
the only enzyme in the pathway that is thought to be 
exclusively committed to TAG synthesis, and thus it is 
considered a key enzyme in this reaction [10,13,16]. 
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Figure 5 Evaluation of triacylglyceride biosynthesis in yeast 
mutants complemented with castor bean DGAT1 or DGT2. The 

H1 122 strain is an are! are2 double mutant, while H1246 is an arel 
ore2 lro1 dgol quadruple mutant. Lipids were separated in a silica- 
gel plate of thin-layer chromatography, and the separation range 
corresponding to Triacylglicerol (TAG) and ADAG (mono-alkyl- 
diacylglicerol), used as lipid standards, is shown after Coomassie- 
Blue R250 staining. 



DGAT activity was first reported in 1956 [65], but it was 
only in the last decade that DGAT genes were cloned 
and molecular tools for studying TG synthesis became 
available [2], Although many studies have been con- 
ducted with these enzymes, little is known about their 
evolutionary origins. Taking advantage of the wide avail- 
ability of data on the genomes of several species, this 
study focused on the characterization of the origin and 
evolutionary history of the DGAT1 and DGAT2 gene 
families in eukaryotes. Our survey of DGAT1 and DGAT2 
proteins shows that they are present in all eukaryotes 
searched, with the exception of yeast (Saccharomyces cere- 
viseae and Candida albicans) and basidiomycetes (Lac- 
caria bicolor, Schizophyllum commune and Agaricus 
bisporus), which lack DGAT1 (Figure 1 and Additional 
Files 1, 2 and 3). This demonstrates that the DGAT1 and 
DGAT2 genes encode essential enzymes for most eukar- 
yotic organisms. In the yeast S. cereviseae, two genes that 
have homology with DGAT1 were identified, ARE1 and 



Eukaryote 
emergence 




Figure 6 Scenarios proposed for the evolution of DGAT1 and 
DGAT2 in eukaryotes. The proposed evolutionary relationships 
between DGAT1 and DGAT2 family members are as follows: DGAT1 
and DGAT2 have evolved asymmetrically since the eukaryotic 
emergence and underwent convergent evolution likely due to the 
importance of TAGs for all eukaryotic organisms. 



ARE2, which encode ACAT-related enzymes. The yeast 
enzymes are 49 percent identical to each other and exhibit 
23 percent identity to human ACAT [66]. Both these 
genes have been suggested to be involved in the synthesis 
of sterol esters. ARE1 is also involved in the synthesis of 
TAG, whereas the ARE2 gene is more specifically involved 
in the synthesis of sterol esters [67]. Insights into the func- 
tions of DGAT1 and DGAT2 in triglyceride metabolism 
have been provided by studies in yeast. With gene knock- 
out or overexpression studies, several groups have demon- 
strated that the yeast DGA1 (DGAT2 homologue), is the 
major DGAT enzyme contributing to triglyceride synthesis 
and storage in yeast [60,68,69]. However, it has been 
demonstrated with genetic complementation that DGAT1 
from Arabidopsis thaliana can restore TAG synthesis to 
the yeast dgal:lrol:arel:are2 quadruple mutant, which is 
devoid of TAG and sterol esters [70]. In the present work, 
using the castor bean genes, we obtained similar results 
with DGAT1, and were also unable to demonstrate 
DGAT2 ability to complement TAG synthesis compared 
to DGAT1, as previously described by Zhang and colla- 
borator's using DG AT2 from Arabidopsis [70] . This indi- 
cates that despite their apparent convergent functional 
evolution, DGAT1 and DGAT2 have different affinities for 
identical substrates, and their effective activity with accu- 
mulation of TAG depends on how much the other 
enzymes from the same pathway consume or provide 
DGAT substrate and precursors. The absence of DGAT1 
genes in yeast and in basidiomycetes indicates that TGAs 
are not essential in these organisms and that other 
enzymes may be able to produce the required amount of 
triacylglycerol. Mutant S. cerevisiae strains devoid of 
DGAT activity and lacking triglycerides are viable and 
grow normally, indicating that triglycerides are not essen- 
tial for the survival of this yeast [60]. However, a recent 
study showed that triglyceride synthesis is essential for via- 
bility of the fission yeast Schizosaccharomyces pombe [71]. 

Our phylogenetic analysis confirmed the existence of 
separate clades for DGAT1 and DGAT2 genes in eukar- 
yotes, indicating that they evolved asymmetrically. In the 
phylogenetic analysis, DGAT1 was more closely related 
and grouped with ACAT genes. ACAT is involved in cho- 
lesterol ester biosynthesis in mammals [28], and homologs 
of these genes were found in fungi, insects and vertebrates 
but not in plants. We verified a conserved motif between 
DGAT1 and ACAT (FYXDWWN; amino acids 392-398 
of A. thaliana DGAT1) in all organisms analyzed that has 
also been reported in other studies [43,72]. This motif has 
been implicated in the binding of fatty acyl-CoA, a com- 
mon substrate for these enzymes [72], Additionally, 
DGAT2 shares homology and groups with acyl-CoA: 
monoacylglycerol acyltransferases [31,73]. Topological stu- 
dies indicated that the majority of DGAT2 is composed of 
a C-terminal cytoplasmic domain. The amino acids 
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HPHG make up residues 161-164 of murine DGAT2, sug- 
gesting that this sequence plays a crucial role in the func- 
tion of the enzyme and may be part of the active site [74]. 
The HPHG-motif is conserved in DGAT2 in yeast, fungi, 
invertebrates and vertebrates and is also present in MGAT 
and AWAT. Plant DGAT2 contains a corresponding 
sequence of conserved amino acids, but the motif is modi- 
fied to EPHSV (108-112 of A thaliana DGAT2), which is 
conserved in all plant species analyzed. Three DGAT2 
homologs were found in the microalgae Chlorella sp. 
(chrDGAT2A, chrDGAT2B and chrDGAT2C). In NJ and 
ML analysis, chrDGAT2C is more related with vertebrate, 
invertebrate, fungi and yeast DGAT2 members and also 
with the MGAT and AWAT proteins. chrDGAT2C has 
an HPHG motif and chrDGAT2A and chrDGAT2B have 
EPHSV motifs, similar to the plant DGAT2. The same 
pattern was found in two unicellular photoautotrophic 
green algae, Ostreococcus tauri and O. lucimarinus, which 
have three DGAT2 genes. One of them shows significant 
homology to mammalian MGAT, whereas the other two 
seem to be more closely related to plant DGAT2 members 
[75]. This suggests that DGAT2 in these organisms may 
have evolved by different duplication events. 

Notably, the DGAT1 and DGAT2 families arose from 
different ancestors during the emergence of eukaryotes, 
and our results suggest that they followed convergent 
evolution in eukaryotes despite having evolved separately 
since the early eukaryotes, as depicted in the model in 
Figure 6. The high sequence similarity and close phyloge- 
netic relationship of the DGAT1 and ACAT genes, 
which is higher than that of DGAT1 and DGAT2, 
strongly suggests that DGAT1 and ACAT emerged from 
a common ancestral sequence, whereas DGAT2 emerged 
from another ancestral sequence. We propose that the 
DGAT1 and DGAT2 genes had a functional convergence 
in eukaryotes, which may be explained by the great 
importance of DGAT activity to the formation of TAGs 
in all eukaryotes. In several species, the formation of 
TAGs plays important roles not only as a storage reserve 
but also in growth and development. In agreement with 
our findings, it has been proposed that the DGAT1 and 
DGAT2 enzymes have different biochemical properties 
and physiological functions. In overexpression studies in 
mice, DGAT2 was identified as a more active enzyme 
than DGAT1, yielding a higher increase in intracellular 
triacylglycerols that accumulated as large, centrally 
located, cytosolic droplets [76]. Mice lacking DGAT1 
(Dgatl" /_ mice) are viable and have modest reductions in 
tissue triacylglycerols [77], whereas mice lacking DGAT2 
(Dgat2" /_ mice) have severe reductions in whole body 
triacylglycerols and die shortly after birth [76]. In A. 
thaliana, mutation of the DGAT1 gene severely affects 
seed development, causing the appearance of wrinkled 
and incompletely filled seeds, aberrant seedling growth, 



reduced seed oil content, an increased DAG/TAG ratio 
and the beta-oxidation of elevated fatty acid (FA) levels 
[13,15]. It has been also demonstrated that DGAT1, but 
not DGAT2, is a multifunctional acyltransferase that cat- 
alyzes the synthesis of diacylglycerol, retinyl esters, and 
waxes in addition to triacylglycerol in in vitro assays [78] . 
DGAT1 cannot compensate for the absence of DGAT2, 
highlighting the functional differences between these 
enzymes. Furthermore, compared to DGAT1, DGAT2 
appears to be more active at lower concentrations of 
acyl-CoA (< 50 [iM) and less active at concentrations 
below 50 mM magnesium [38]. The two enzymes may 
also differ in terms of substrate specificities. Substrate 
specificity assays using a mammalian DGAT2 expressed 
in insect cells showed that oleoyl-CoA was utilized pre- 
ferentially, followed by palmitoyl-CoA, while activity 
levels with linoleoyl- and arachidonyl-CoA were similar 
[38]. However, in enzymatic assays using soluble and 
non-soluble samples of the fungal DGAT2 from M. 
ramanniana, no significant preference was detected 
between 12:0-CoA and 18:l-CoA as an acyl donor sub- 
strate [30]. In plant species, although both DGAT1 and 
DGAT2 significantly contribute to TAG formation, the 
relative contribution of each enzyme seems to be species 
dependent and may also differ in different tissues within 
the same plant species [40]. Moreover, it was reported 
that DGAT1 and DGAT2 are located in different areas of 
the endoplasmic reticulum (ER) [40] . The fatty acid pro- 
file of TAG may be determined by distinct substrate spe- 
cificities of the DGAT enzymes in different species and 
in plants containing unusual fatty acids. DGAT2 may 
have a major role in channeling unusual fatty acids into 
seed storage oils [15,40]. In a recent study with Ricinus 
comunis, it was demonstrated that DGAT2 presents 
higher mRNA accumulation than DGAT1 during castor 
seed development [79]. It was also demonstrated in R. 
comunis that only RcDGAT2 can increase hydroxy fatty 
acid levels in transgenic Arabidopsis, while RcDGATl 
cannot [80]. DGAT2 is more likely to play a major role in 
seed TAG biosynthesis than DGAT1 in this plant species 
[44]. In addition, a study of Vernolic acid in Vernonia 
galamensis concluded that both DGAT1 and DGAT2 
increase epoxy fatty acid accumulations, but DGAT2 had 
a greater effect [81]. All these findings corroborate with 
our results that show that DGAT1 and DGAT2, despite 
being evolutionarily divergent enzymes, play important 
roles in TAG metabolism in all eukaryotic organisms. 
Therefore, we conclude that they are very ancient 
enzymes that arose quite early in eukaryotic evolution. 

The same patterns were observed in structural and 
transmembrane analysis. The topology models based on 
prediction algorithms suggested that DGAT1 and ACAT 
had eight transmembrane domains and 15-17 exons, 
whereas DGAT2 had 1 or 2 transmembrane domains and 
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8 exons, similar to MGAT and AWAT. This reinforced 
the idea that despite their ability to catalyze similar reac- 
tions, DGAT1 and DGAT2 belong to different gene 
families that share neither DNA nor protein sequence 
similarity [2]. In addition, topological studies show that 
DGAT1 has its N- terminus oriented towards the cytosol 
and C-terminal region, which accounts for approximately 
50% of the protein, and is present in the ER lumen [82]. 
Conversely, DGAT2 is an integral membrane protein 
with both the N and C termini oriented toward the cyto- 
sol [74]. In agreement with our results, a comparative 
genomics and proteomics study of vertebrate DGATs 
revealed that the predicted secondary and transmem- 
brane (TrM) structures of DGAT1 and DGAT2 were dis- 
tinct [50]. All of the DGAT2 subunits had one or two 
predicted TrM(s), whereas DGAT1 was predicted to have 
nine regions with high probability for TrM [50]. In con- 
trast, a recent study about the membrane topology of 
murine DGAT1 indicates that DGAT1 has three trans- 
membrane domains [82]. For murine DGAT2, two trans- 
membrane domain were identified [74] . Studies were also 
carried out in plants and have shown that Helianthus 
annuus [48], Echium [83] and Brassica juncea [43] have 
nine transmembrane domains in DGAT1, while DGAT2 
from Brassica juncea have two transmembrane domains 
[43]. 

The Ka/Ks ratios for vertebrate DGAT1, DGAT2, 
ACAT, MGAT and AWAT and for plant DGAT1 and 
DGAT2 showed that distinct selective pressures have 
been operating on the DGAT1 and DGAT2 genes. The 
comparisons between species groups within of each 
DAGT1 and DAGT2 family revealed high /Cy/^ratios, 
while low K a /K s ratios were observed between the 
DGAT1 and DGAT2 families. This indicates a purifying 
selection within each family with a strong positive selec- 
tion between them. In the DGAT1 and ACAT genes, in 
addition to a high level of similarity between them, a 
positive selection has taken place at multiple locations. 

Conclusions 

Overall, our data show that the DGAT1 and DGAT2 
proteins are present in most eukaryotic organisms as 
different gene families. Phylogenetic and evolutionary 
analyses revealed that DGAT1 and DGAT2 evolved 
separately with functional convergence since the emer- 
gence of eukaryotes and are widely divergent in molecu- 
lar and structural aspects. 

Additional material 



Additional file 1: The DGAT1 , DGAT2, ACAT, MGAT, AWAT and 
other enzymes used for the analyses. The table shows the species, 
gene names and accession numbers of the sequences used in the 
analyses. 



Additional file 2: Phylogenetic tree of DGAT1 and DGAT2 gene 
families reconstructed by the Neighbor joining (NJ) method. A total 
of 151 protein sequences from eukaryotic species and 344 sites were 
included in the analyses. The bootstrap values are labeled above the 
branches. Only values > 50% are presented. Numbers within brackets 
correspond to the number of species within each group. The enzymes 
are represented by different colors in the phylogenetic tree: DGAT1 (red), 
DGAT2 (green), ACATs (purple), MGAT (dark blue) and AWAT (light blue). 
Taxa terminologies are abbreviated using the first letter of the genus and 
two letters of the species name: Alternaha brassicicola (Abr), Arabidopsis 
thaliana (Ath), Dictyostelium discoideum (Ddi), Euonymus alatus (Eal), 
Homo sapiens (Hsa), Mortierella ramanniana (Mra), Thalassiosira 
pseudonana (Tps), Xenopus tropicalis (Xtr). 

Additional file 3: Phylogenetic tree of the DGAT1 and DGAT2 gene 
families reconstructed by the Maximum likelihood (ML) method. A 

total of 151 protein sequences from eukaryotic species and 344 sites 
were included in the analyses. The bootstraps values are labeled above 
the branches. Only values > 50% are presented. Numbers within brackets 
correspond to the number of species within each group. The enzymes 
are represented by different colors in the phylogenetic tree: DGAT1 (red), 
DGAT2 (green), ACATs (purple), MGAT (dark blue) and AWAT (light blue). 
Taxa terminologies are abbreviated using the first letter of the genus and 
two letters of the species name: Alternaria brassicicola (Abr), Arabidopsis 
thaliana (Ath), Dictyostelium discoideum (Ddi), Euonymus alatus (Eal), 
Homo sapiens (Hsa), Mortierella ramanniana (Mra), Thalassiosira 
pseudonana (Tps), Xenopus tropicalis (Xtr). 

Additional file 4: Multiple sequence alignment of predicted amino 
acid sequences of DGAT1 proteins. DGAT1 sequences from A. thaliana 
(Ath), M. truncatula (Mtr), R. comunis (Rco), 0. sativa (Osa), Z mays (Zma) 
and M. musculus (Mmu), and ACAT proteins from M. musculus and 5. 
cerevisieae (See) were aligned. Identical residues are shaded black, and 
similar residues are shaded gray. The DAG/phorbol ester binding 
signature motif is underlined with a dotted line; the triangle (▼) shows 
the conserved phenylalanine. The conserved motif between DGAT1 and 
ACATs is underlined. 

Additional file 5: Multiple sequence alignment of deduced amino 
acid sequences of DGAT2 proteins. DGAT2 sequences from A. thaliana 
(Ath), M. truncatula (Mtr), R. comunis (Rco), 0. sativa (Osa), Z mays (Zma), 
M. musculus (Mmu) and 5. cerevisieae (See) were aligned with MGAT and 
AWAT proteins from M. musculus. Identical residues are shaded in black, 
and similar residues are shaded in gray. 

Additional file 6: Estimation of Ka/Ks rates on vertebrate sequences. 

Substitution rates between full-length cDNA sequences of vertebrate 
DGAT1, DGAT2, ACAT1, ACAT2, MGAT and AWAT. Comparisons were 
performed between ACATs and DGAT1, MGATs and DGAT2 and between 
AWAT and DGAT2. 

Additional file 7: Estimation of Ka/Ks rates on plant sequences. 

Substitution rates between full-length cDNA sequences of plant DGAT1 
and DGAT2. Comparisons were performed for DGAT1 and DGAT2 
between monocotyledon and eudicotyledon, between oilseed and non- 
oilseed plants and between DGAT1 and DGAT2 from all plant species. 



Acknowledgements 

This work was supported by CNPq (Conselho Nacional de Desenvolvimento 
Cientifico e Tecnologico) grant number 559636/2009-1, CAPES (Coordenacao 
de Aperfeicoamento de Pessoal de Nivel Superior), FAPERGS (Fundacao de 
Amparo a Pesquisa do Estado do Rio Grande do Sul), FINEP (Financiadora de 
Projetos) and MCT (Ministerio de Ciencia e Tecnologia). A. Turchetto-Zolet 
and F. Maraschin have PNPD/CAPES fellowships, G. Loss has a Ph.D. 
fellowship from CAPES, and A. Cagliari received a Ph.D. fellowship from 
CNPq. M. Margis-Pinheiro and R. Margis are recipients of CNPq research 
fellowships numbers 308708/2006-7 and 303967/2008-0, respectively. 

Author details 

1 Programa de P6s-Graduacao em Genetica e Biologia Molecular, 
Departamento de Genetica, Universidade Federal do Rio Grande do Sul, 
Brazil. 2 Centro de Biotecnologia, Universidade Federal do Rio Grande do Sul, 



Turchetto-Zolet et al. BMC Evolutionary Biology 201 1, 1 1:263 
http://www.biomedcentral.eom/1471-2148/1 1/263 



Page 12 of 14 



Brazil. 3 Departamento de Biofisica, Universidade Federal do Rio Grande do 
Sul, Brazil. 

Authors' contributions 

ACTZ retrieved the sequences used in the study, made the sequence 
alignments and performed the phylogenetics, structural and evolution 
analyses. FSM contributed to cloning and yeast complementation. RM and 
CMBA performed the lipid analysis. MMP and RM conceived the study. RM, 
MMP, GLM, and AC contributed to the analysis interpretation. All authors 
contributed to the writing of the manuscript and read and approved the 
final manuscript. 

Received: 31 January 201 1 Accepted: 20 September 201 1 
Published: 20 September 201 1 

References 

1. Athenstaedt K, Daum G: The life cycle of neutral lipids: synthesis, storage 

and degradation. Cell Mol Life Sci 2006, 63(1 2):1 355-1 369. 

2. Yen CLE, Stone SJ, Koliwad S, Harris C, Farese RV: DGAT enzymes and 
triacylglycerol biosynthesis. J Lipid Res 2008, 49(1 1):2283-2301. 

3. Murphy DJ: Plant Lipids: Biology, Utilization and Manipulation Oxford: 
Blackwell Publishing; 2005. 

4. Durrett TP, Benning C, Ohlrogge J: Plant triacylglycerols as feedstocks for 
the production of biofuels. Plant J 2008, 54(4):593-607. 

5. Dyer JM, Stymne S, Green AG, Carlsson AS: High-value oils from plants. 
Plant J 2008, 54(4):640-655. 

6. Voelker T, Kinney AT: Variations in the biosynthesis of seed-storage lipids. 
Annual Review of Plant Physiology and Plant Molecular Biology 2001, 
52:335-361. 

7. Huang AHC: Oil bodies and oleosins in seeds. Annual Review of Plant 
Physiology and Plant Molecular Biology 1992, 43:1 77-200. 

8. Kennedy EP, Weiss SB: Function of cytidine coenzymes in the 
biosynthesis of phospholipides. J Biol Chem 1956, 222(1 ):1 93-21 3. 

9. Ohlrogge J, Browse J: Lipid biosynthesis. Plant Cell 1995, 7(7):957-970. 

10. Ichihara K, Takahashi T, Fujii S: Diacylglycerol Acyltransferase in Maturing 
Safflower Seeds-Its Influences on the Fatty-Acid Composition of 
Triacylglycerol and on the Rate of Triacylglycerol Synthesis. Biochimica Et 
Biophysica Acta 1988, 958(1 ):1 25-1 29. 

11. Perry HJ, Harwood JL: Changes in the lipid content of developing seeds 
of Brassica napus. Phytochemistry 1993, 32(6):1 41 1-1415. 

12. Lardizabal K, Effertz R, Levering C, Mai J, Pedroso MC, Jury T, Aasen E, 
Gruys K, Bennett K: Expression of Umbelopsis ramanniana DGAT2A in 
seed increases oil in soybean. Plant Physiol 2008, 148(1 ):89-96. 

13. Lung SC, Weselake RJ: Diacylglycerol acyltransferase: A key mediator of 
plant triacylglycerol synthesis. Lipids 2006, 41 (12):1 073-1 088. 

14. Slabas AR, Simon JW, Brown AP: Biosynthesis and regulation of fatty acids 
and triglycerides in oil seed rape. Current status and future trends. Eur J 
Lipid Sci Tech 2001, 103(7):455-466. 

15. Xu JY, Francis T, Mietkiewska E, Giblin EM, Barton DL, Zhang Y, Zhang M, 
Taylor DC: Cloning and characterization of an acyl-CoA-dependent 
diacylglycerol acyltransferase 1 (DGAT1) gene from Tropaeolum majus, 
and a study of the functional motifs of the DGAT protein using site- 
directed mutagenesis to modify enzyme activity and oil content. Plant 
Biotechnol J 2008, 6(8):799-818. 

16. Jako C, Kumar A, Wei YD, Zou JT, Barton DL, Giblin EM, Covello PS, 
Taylor DC: Seed-specific over-expression of an Arabidopsis cDNA 
encoding a diacylglycerol acyltransferase enhances seed oil content and 
seed weight. Plant Physiol 2001, 1 26(2):861 -874. 

17. Settlage SB, Kwanyuen P, Wilson RF: Relation between diacylglycerol 
acyltransferase activity and oil concentration in soybean. J Am Oil Chem 
Soc 1998, 75(7):775-781. 

18. Lock YY, Snyder CL, Zhu WM, Siloto RMP, Weselake RJ, Shah S: Antisense 
suppression of type 1 diacylglycerol acyltransferase adversely affects 
plant development in Brassica napus. Physiologia Plantarum 2009, 
137(1):61-71. 

19. Zheng P, Allen WB, Roesler K, Williams ME, Zhang S, Li J, Glassman K, 
Ranch J, Nubel D, Solawetz W, et al: A phenylalanine in DGAT is a key 
determinant of oil content and composition in maize. Nat Genet 2008, 
40(3):367-372. 

20. Dahlqvist A, Stahl U, Lenman M, Banas A, Lee M, Sandager L, Ronne H, 
Stymne H: Phospholipid: diacylglycerol acyltransferase: An enzyme that 



catalyzes the acyl-CoA-independent formation of triacylglycerol in yeast 
and plants. P Natl Acad Sci USA 2000, 97(12):6487-6492. 

21. Stahl U, Carlsson AS, Lenman M, Dahlqvist A, Huang BQ, Banas W, Banas A, 
Stymne S: Cloning and functional characterization of a Phospholipid: 
Diacylglycerol acyltransferase from Arabidopsis. Plant Physiol 2004, 
135(3):1 324-1 335. 

22. Banas A, Dahlqvist A, Stahl U, Lenman M, Stymne S: The involvement of 
phospholipid: diacylglycerol acyltransferases in triacylglycerol 
production. Biochem Soc T 2000, 28:703-705. 

23. Kalscheuer R, Steinbuchel A: A novel bifunctional wax ester synthase/acyl- 
CoA: diacylglycerol acyltransferase mediates wax ester and 
triacylglycerol biosynthesis in Acinetobacter calcoaceticus ADP1. J Biol 
Chem 2003, 278(1 0):8075-8082. 

24. King A, Nam JW, Han JX, Hilliard J, Jaworski JG: Cuticular wax biosynthesis 
in Petunia petals: cloning and characterization of an alcohol- 

acy transferase that synthesizes wax-esters. Planta 2007, 226(2):381 -394. 

25. Li F, Wu X, Lam P, Bird D, Zheng H, Samuels L, Jetter R, Kunst L: 
Identification of the wax ester synthase/acyl-coenzyme A: Diacylglycerol 
acyltransferase WSD1 required for stem wax ester biosynthesis in 
Arabidopsis. Plant Physiol 2008, 1 48(1 ):97-1 07. 

26. Saha S, Enugutti B, Rajakumari S, Rajasekharan R: Cytosolic triacylglycerol 
biosynthetic pathway in oilseeds. Molecular cloning and expression of 
peanut cytosolic diacylglycerol acyltransferase. Plant Physiol 2006, 

1 41 (4):1 533-1 543. 

27. Durrett TP, McClosky DD, Tumaney AW, Elzinga DA, Ohlrogge J, Pollard M: 
A distinct DGAT with sn-3 acetyltransferase activity that synthesizes 
unusual, reduced-viscosity oils in Euonymus and transgenic seeds. P Natl 
Acad Sci USA 2010, 107(20):9464-9469. 

28. Cases S, Smith SJ, Zheng YW, Myers HM, Lear SR, Sande E, Novak S, 
Collins C, Welch CB, Lusis AJ, et al: Identification of a gene encoding an 
acyl CoA:diacylglycerol acyltransferase, a key enzyme in triacylglycerol 
synthesis. P Natl Acad Sci USA 1998, 95(22):1 301 8-1 3023. 

29. Hofmann K: A superfamily of membrane-bound O-acyltransferases with 
implications for Wnt signaling. Trends Biochem Sci 2000, 25(3):1 11-112. 

30. Lardizabal KD, Mai JT, Wagner NW, Wyrick A, Voelker T, Hawkins DJ: DGAT2 
is a new diacylglycerol acyltransferase gene family-Purification, cloning, 
and expression in insect cells of two polypeptides from Mortierella 
ramanniana with diacylglycerol acyltransferase activity. J Biol Chem 2001, 
276(42):38862-38869. 

31. Yen CLE, Stone SJ, Cases S, Zhou P, Farese RV: Identification of a gene 
encoding MGAT1, a monoacylglycerol acyltransferase. P Natl Acad Sci 
USA 2002, 99(1 3):85 12-85 17. 

32. Cao JS, Lockwood J, Burn P, Shi YG: Cloning and functional 
characterization of a mouse intestinal Acyl-CoA: monoacylglycerol 
acyltransferase, MGAT2. J Biol Chem 2003, 278(1 6):1 3860-1 3866. 

33. Yen CLE, Farese RV: MGAT2, a monoacylglycerol acyltransferase 
expressed in the small intestine. J Biol Chem 2003, 278(20):1 8532-1 8537. 

34. Cheng D, Nelson TC, Chen J, Walker SG, Wardwell-Swanson J, Meegalla R, 
Taub R, Billheimer JT, Ramaker M, Feder JN: Identification of acyl 
coenzyme A: Monoacylglycerol acyltransferase 3, an intestinal specific 
enzyme implicated in dietary fat absorption. J Biol Chem 2003, 

278(1 6):1 361 1-13614. 

35. Cheng JB, Russell DW: Mammalian wax biosynthesis-l. Identification of 
two fatty acyl-coenzyme A reductases with different substrate 
specificities and tissue distributions. J Biol Chem 2004, 
279(36):37789-37797. 

36. Turkish AR, Henneberry AL, Cromley D, Padamsee M, Oelkers P, Bazzi H, 
Christiano AM, Billheimer JT, Sturley SL: Identification of two novel human 
acyl-CoA wax alcohol acyltransferases: members of the diacylglycerol 
acyltransferase 2 (DGAT2) gene superfamily. J Biol Chem 2005, 

280(1 5):1 4755-1 4764. 

37. Cases S, Smith SJ, Zheng Y, Myers HM, Sande ER, Novak S, Lear SR, 
Erickson SK, Farese RV: Cloning and expression of a candidate gene for 
diacylglycerol acyltransferase. Faseb J 1998, 1 2(5):A81 4-A81 4. 

38. Cases S, Stone SJ, Zhou P, Yen E, Tow B, Lardizabal KD, Voelker T, Farese RV: 
Cloning of DGAT2, a second mammalian diacylglycerol acyltransferase, 
and related family members. J Biol Chem 2001, 276(42):38870-38876. 

39. Chen GQ, Turner C, He X, Nguyen T, McKeon TA, Laudencia- 
Chingcuanco D: Expression profiles of genes involved in fatty acid and 
triacylglycerol synthesis in castor bean {Ricinus communis L). Lipids 2007, 
42(3):263-274. 



Turchetto-Zolet et al. BMC Evolutionary Biology 201 1, 1 1:263 
http://www.biomedcentral.eom/1471-2148/1 1/263 



40. Shockey JM, Gidda SK, Chapital DC, Kuan JC, Dhanoa PK, Bland JM, 
Rothstein SJ, Mullen RT, Dyer JM: Tung tree DGAT1 and DGAT2 have 
nonredundant functions in triacylglycerol biosynthesis and are localized 
to different subdomains of the endoplasmic reticulum. Plant Cell 2006, 

18(9):2294-2313. 

41. Banilas G, Karampelias M, Makariti I, Kourti A, Hatzopoulos P: The olive 
DGAT2 gene is developmentally regulated and shares overlapping but 
distinct expression patterns with DGAT1 . J Exp Bot 201 1. 

42. He XH, Chen GQ, Lin JT, McKeon TA: Regulation of diacylglycerol 
acyltransferase in developing seeds of castor. Lipids 2004, 39(9):865-871. 

43. Ilaiyaraja N, Rajarani AP, Santha IM: Cloning and characterization of 
diacylglycerol acyltransferase (DGAT) cDNA sequence from Brossico 
juncea cv Pusa Bold. Indian Journal of Biochemistry & Biophysics 2008, 
45(1):30-36. 

44. Kroon JTM, Wei WX, Simon WJ, Slabas AR: Identification and functional 
expression of a type 2 acyl-CoA: diacylglycerol acyltransferase (DGAT2) 
in developing castor bean seeds which has high homology to the major 
triglyceride biosynthetic enzyme of fungi and animals. Phytochemistry 

2006, 67(23):2541-2549. 

45. Milcamps A, Tumaney AW, Paddock T, Pan DA, Ohlrogge J, Pollard M: 
Isolation of a gene encoding a 1,2-diacylglycerol-sn-acetyl-CoA 
acetyltransferase from developing seeds of Euonymus alatus. J Biol Chem 
2005, 280(7):5370-5377. 

46. Nykiforuk CL, Furukawa-Stoffer TL, Huff PW, Sarna M, Laroche A, 
Moloney MM, Weselake RJ: Characterization of cDNAs encoding 
diacylglycerol acyltransferase from cultures of Brossica napus and 
sucrose-mediated induction of enzyme biosynthesis. BBA-Mol Cell Biol L 
2002, 1580(2-3):95-109. 

47. Routaboul JM, Benning C, Bechtold N, Caboche M, Lepiniec L: The TAG1 
locus of Arabidopsis encodes for a diacylglycerol acyltransferase. Plant 
Physiology and Biochemistry 1999, 37(1 1):83 1-840. 

48. Sun L, Ouyang C, Kou SL, Wang SH, Yao YY, Peng T, Xu Y, Tang L, Chen F: 
Cloning and Characterization of a cDNA Encoding Type 1 Diacylglycerol 
Acyltransferase from Sunflower {Helianthus annuus L). Z Naturforsch C 
2011, 66(1-2):63-72. 

49. Wang HW, Zhang JS, Gai JY, Chen SY: Cloning and comparative analysis 
of the gene encoding diacylglycerol acyltransferase from wild type and 
cultivated soybean. Theor App Genet 2006, 1 12(6):1 086-1 097. 

50. Holmes RS: Comparative genomics and proteomics of vertebrate 
diacylglycerol acyltransferase (DGAT), acyl CoA wax alcohol 
acyltransferase (AWAT) and monoacylglycerol acyltransferase (MGAT). 
Comp Biochem Physiol Part D Genomics Proteomics 2010, 5(1 ):45-54. 

51. Andrianov V, Borisjuk N, Pogrebnyak N, Brinker A, Dixon J, Spitsin S, Flynn J, 
Matyszczuk P, Andryszak K, Laurelli M, et al: Tobacco as a production 
platform for biofuel: overexpression of Arabidopsis DGAT and LEC2 
genes increases accumulation and shifts the composition of lipids in 
green biomass. Plant Biotechnol J 2010, 8(3):277-287. 

52. Klaus D, Ohlrogge JB, Neuhaus HE, Dormann P: Increased fatty acid 
production in potato by engineering of acetyl-CoA carboxylase. Planta 
2004, 219(3)389-396. 

53. Thompson JD, Higgins DG, Gibson TJ: Clustal-W-lmproving the sensitivity 
of progressive multiple sequence alignment through sequence 
weighting, position-specific gap penalties and weight matrix choice. 
Nucleic Acids Res 1 994, 22(22):4673-4680. 

54. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary 
genetics analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 
24(8):1 596-1 599. 

55. Stoveken T, Kalscheuer R, Malkus U, Reichelt R, Steinbuchel A: The wax 
ester synthase/acyl coenzyme A: diacylglycerol acyltransferase from 
Acinetobacter sp strain ADP1: Characterization of a novel type of 
acyltransferase. J Bacteriol 2005, 187(4):1 369-1 376. 

56. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate 
large phylogenies by maximum likelihood. Syst Biol 2003, 52(5):696-704. 

57. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference 
under mixed models. Bioinformatics 2003, 1 9(1 2):1 572-1 574. 

58. Solovyev W: Statistical approaches in eukaryotic gene prediction. In 
Handbook of statistical genetics. Edited by: Balding D. New York: John Wiley 
2001:83-127. 

59. Librado P, Rozas J: dnaSP v5: a software for comprehensive analysis os 
DNA polimorphism data. Bioinformatics 2009, 25(1451-1452). 



Page 13 of 14 



60. Sandager L, Gustavsson MH, Stahl U, Dahlqvist A, Wiberg E, Banas A, 
Lenman M, Ronne H, Stymne S: Storage lipid synthesis is non-essential in 
yeast. J Biol Chem 2002, 277(8):6478-6482. 

61. Mangold HK, Malins DC: Fractions of fats, oils, and waxes on thin layer of 
silic acid. J Am Oil Chem Soc 1960, 37(8)383-385. 

62. Nakamura K, Handa S: Coomassie Brilliant Blue staining of lipids on thin- 
layer plates. Anal Biochem 1984, 142(2):406-410. 

63. Nei M, Gojobori T: Simple methods for estimating the numbers of 
synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol 
1986, 3(5):41 8-426. 

64. Yang Z, Bielawski JP: Statistical methods for detecting molecular 
adaptation. Trends Ecol Evol 2000, 1 5(1 2):496-503. 

65. Weiss SB, Kennedy EP, Kiyasu JY: The enzymatic synthesis of triglycerides. 
J Biol Chem 1 960, 235:40-44. 

66. Yang HY, Bard M, Bruner DA, Gleeson A, Deckelbaum RJ, Aljinovic G, 

Pohl TM, Rothstein R, Sturley SL: Sterol esterification in yeast: A two-gene 
process. Science 1996, 272(5266):1 353-1 356. 

67. Sandager L, Dahlqvist A, Banas A, Stahl U, Lenman M, Gustavsson M, 
Stymne S: An acyl-CoA: cholesterol acyltransferase (ACAT)-related gene is 
involved in the accumulation of triacylglycerols in Saccharomyces 
cerevisiae. Biochem Soc T 2000, 28:700-702. 

68. Oelkers P, Cromley D, Padamsee M, Billheimer JT, Sturley SL: The DGA1 
gene determines a second triglyceride synthetic pathway in yeast. J Biol 
Chem 2002, 277(1 1):8877-8881. 

69. Sorger D, Daum G: Synthesis of triacylglycerols by the acyl-coenzyme A: 
Diacyl-glycerol acyltransferase Dgalp in lipid particles of the yeast 
Saccharomyces cerevisiae. J Bacteriol 2002, 184(2):5 19-524. 

70. Zhang M, Fan JL, Taylor DC, Ohlrogge JB: DGAT1 and PDAT1 
Acyltransferases Have Overlapping Functions in Arabidopsis 
Triacylglycerol Biosynthesis and Are Essential for Normal Pollen and 
Seed Development. Plant Cell 2009, 21(12)3885-3901. 

71. Zhang Q, Chien HK, Low C, Zhang SC, Yang HY: Synthesis of 
triacylglycerols is required for the fission yeast Schizosaccharomyces 
pombe to enter stationary phase. Yeast 2003, 20:S297-S297. 

72. Guo ZM, Cromley D, Billheimer JT, Sturley SL: Identification of potential 
substrate-binding sites in yeast and human acyl-CoA sterol 
acyltransferases by mutagenesis of conserved sequences. J Lipid Rese 
2001, 42(8):1 282-1 291. 

73. Lockwood JF, Cao JS, Burn P, Shi YG: Human intestinal monoacylglycerol 
acyltransferase: differential features in tissue expression and activity. Am 
J Physiol-Endoc M 2003, 285(5):E927-E937. 

74. Stone SJ, Levin MC, Farese RV: Membrane topology and identification of 
key functional amino acid residues of murine Acyl-CoA: Diacylglycerol 
acyltransferase-2. J Biol Chem 2006, 281 (52):40273-40282. 

75. Wagner M, Hoppe K, Czabany T, Heilmann M, Daum G, Feussner I, Fulda M: 
Identification and characterization of an acyl-CoA:diacylglycerol 
acyltransferase 2 (DGAT2) gene from the microalga O. tauri. Plant Physiol 
Biochem 2010, 48(6):407-416. 

76. Stone SJ, Myers HM, Watkins SM, Brown BE, Feingold KR, Elias PM, 
Farese RV: Lipopenia and skin barrier abnormalities in DGAT2-deficient 
mice. J Biol Chem 2004, 279(1 2):1 1 767-1 1 776. 

77. Smith SJ, Cases S, Jensen DR, Chen HC, Sande E, Tow B, Sanan DA, Raber J, 
Eckel RH, Farese RV Jr: Obesity resistance and multiple mechanisms of 
triglyceride synthesis in mice lacking Dgat. Nat Genet 2000, 25(1 ):87-90. 

78. Yen CLE, Monetti M, Burri BJ, Farese RV: The triacylglycerol synthesis 
enzyme DGAT1 also catalyzes the synthesis of diacylglycerols, waxes, 
and retinyl esters. J Lipid Res 2005, 46(7):1 502-1 51 1. 

79. Cagliari A, Margis-Pinheiroa M, Loss G, Mastroberti AA, Mariath JEA, 
Margis M: Identification and expression analysis of castor bean {Ricinus 
communis) genes encoding enzymes from the triacylglycerol 
biosynthesis pathway. Plant Sci 2010. 

80. Burgal J, Shockey J, Lu CF, Dyer J, Larson T, Graham I, Browse J: Metabolic 
engineering of hydroxy fatty acid production in plants: RcDGAT2 drives 
dramatic increases in ricinoleate levels in seed oil. Plant Biotechnol J 
2008, 6(8):81 9-831. 

81. Li RZ, Yu KS, Hatanaka T, Hildebrand DF: Vernonia DGATs increase 
accumulation of epoxy fatty acids in oil. Plant Biotechnol J 2010, 
8(2):184-195. 

82. McFie PJ, Stone SL, Banman SL, Stone SJ: Acyl CoA:Diacylglycerol 
acyltransferase-1 (DGAT1): topological orientation, identification of a 



Turchetto-Zolet et al. BMC Evolutionary Biology 201 1, 1 1:263 
http://www.biomedcentral.eom/1471-2148/1 1/263 



Page 14 of 14 



putative active site histidine and the role of the N-terminus in dimer/ 
tetramer formation. J Biol Chem 2010. 
83. Manas-Fernandez A, Vilches-Ferron M, Garrido-Cardenas JA, Belarbi EH, 
Alonso DL, Garcia-Maroto F: Cloning and Molecular Characterization of 
the Acyl-CoA:Diacylglycerol Acyltransferase 1 (DGAT1) Gene from 
Echium. Lipids 2009, 44(6):555-568. 



doi:1 0.1 1 86/1 471 -21 48-1 1 -263 

Cite this article as: Turchetto-Zolet et al:. Evolutionary view of acyl-CoA 
diacylglycerol acyltransferase (DGAT), a key enzyme in neutral lipid 
biosynthesis. BMC Evolutionary Biology 201 1 11 :263. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 

Submit your manuscript at (^\ RioMM i r pnt ^\ 

www.biomedcentral.com/submit Blomea ^ enirai 



